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ABSTRACT 


This  paper  is  concerned  with  the  preparation  of  a  "mathematical 
model"  Involving  collisiona  between  two  bodiea  and  the  subaequient 
analog  computer  program.  It  conaidera  a  flctltloua  damped  aprlng 
aa  a  coefficient  of  reatitution  "generator."  Formulae  are  derived 
for  spring  atiffneaa  and  the  viacoua  damping  in  both  one  and  two 
degreea  of  freedom. 
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SYMBOLISM 

Xi 

displacement  of  1  th  mass 

X  -  XJ-X2 

relative  displacement  of  the  masses 

X 

computer  voltage  analogous  to  x 

Vi 

velocity  of  1  th  mass  just  before  contact 

n 

velocity  of  i  th  mass  just  after  separation 

*i 

velocity  of  1  th  mass  as  a  function  of  time 

X  -  (ki-X2) 

relative  velocity  of  the  masses 

Xi 

acceleration  of  i  th  mass 

t 

time 

tl 

time  when  (X]^-*X2)  is  at  a  maximum 

time  when  (X]^-X2)  ■  0  the  second  time 

machine  time 

M 

mass 

Ml 

mass  of  first  body 

M2 

mass  of  second  body 

R 

damping  factor  (#sec/in) 

K 

s  pr ing  s  t i f  f nes  s  (#/ in) 

a 

coefficient  of  restitution 

q 

naas  ratio 

Primes  denote  velocity  and  displacement  immediately  after  separation, 
i.e.,  relative  displacement  is  zero. 
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Representation  of  the  Coefficient  of  Restitution  In  a  Dynamic 
System  with  Analog  Computer  Applications. 

OBJECTIVES 


1.  To  determine  a  system  of  differential  equations  whereby  the 

<  energy  losses  due  to  the  coefficient  of  restitution  may  be  represented 

as  a  continuous  function  of  time. 

2.  To  formulate  an  Analog  Computer  program  which  simulates 
the  coefficient  of  restitution. 

3.  To  evaluate  the  parameters  representing  the  collisions  be¬ 
tween  barrel  extension  and  buffer,  and  barrel  extension  and  receiver 
in  the  M-73  machine  gun. 

CONCLUSION 

It  has  been  analytically  determined  that  a  fictitious  damped 
spring  produces  energy  losses  analogous  to  those  produced  by  the 
coefficient  of  restitution. 

Formulae  have  been  developed  to  permit  the  evaluation  of  the 
spring  stiffness  and  viscous  damping  factors  as  functions  of  the 
coefficient  of  restitution  and  the  Impact  time.  These  formulae  have 
been  applied  to  the  two  major  collisions  which  occur  In  the  M-73 
automatic  weapon. 

The  Analog  Computer  was  then  used  to  simulate  these  collisions 
and  the  resulting  solution  was  compared  with  test  data.  In  com¬ 
paring  the  results.  It  was  found  that  the  differences  were  those 
due  to  human  error  In  interpretation  of  the  respective  recordings. 
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1.  INTRODUCTION 


a.  When  two  bodies  collide,  a  transfer  of  energy  takes  place 
through  the  active  and  reactive  forces  of  Impact.  These  forces  are 
dependent  on  the  shape  and  elastic  properties  of  the  bodies.  As  a 
rule,  the  energy  before  Impact  Is  greater  than  the  energy  after 
Impact.  The  loss  of  energy  In  the  form  of  heat  and  deformation  may 
by  expressed  by  the  coefficient  of  restitution. 

b.  If  the  time  of  Impact  and  the  elastic  deformation  are  not 
considered,  the  relationships^  between  the  Initial  velocities 
(Immediately  before  contact)  and  the  final  velocities  (Immediately 
after  separation)  are  defined  as  follows: 

/  C'+a)+ 

M.V,  Wz)V;.  A 

^  M,  W*. 

(.V,  - 

where  "a"  Is  defined  as  the  coefficient  of  restitution. 

c.  When  solid  stops  or  buffers  must  be  represented  by  means 
of  dynamic  equations,  serious  errors  may  be  Introduced  Into  the  equa¬ 
tions  of  motion  If  the  collisions  are  assumed  to  be  completely  plastic 
or  elastic.  Therefore,  It  Is  necessary  to  Introduce  realistic  energy 
losses. 


V, 


d.  Experimental  results  show  that  when  two  bodies  collide,  a 
deformation  of  the  bodies  occurs  during  contact  time.  The  bodies  may 
resume  their  original  shape  after  contact.  Consequently,  the  equations 
of  motion  should  be  In  such  a  form  as  to  permit  slight  displacements. 
This  Is  easily  accomplished  by  Introducing  a  fictitious  spring,  having 
a  high  spring  rate.  Into  the  system.  Since  It  Is  necessary  for  the 
spring  to  return  less  energy  Into  the  system  than  It  absorbs  --  viscous 
damping  Is  employed  to  produce  the  required  energy  losses. 


1 


Superscripts  denote  reference  number. 
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e.  Paragraph  (d)  describes  Che  applied  logic  used  Co  develop 
Che  dynamic  equaclons  of  moClon  which  consider  Che  coefflclenC  of 
resClCuClon.  IC  Is  now  possible  Co  conscrucc  an  analog  compuCer 
paCch  diagram  Chrough  vdilch  Che  sysCem  may  be  evaluaCed.  SubsequenC 
secclons  of  Chls  reporC  will: 

(1)  Prove  Che  valldlCy  of  Che  applied  logic; 

(2)  Show  conslsCency  beCween  dynamic  and  sCatlc 
equaclons; 

(3)  Demons Crace  Che  applies Cion  of  Che  dynamic 
equaclons  of  moclon  Co  analogous  syscems  of 
Che  M-73  machine  gun. 

f.  In  Che  course  of  Chls  analysis  1C  Is  assumed  chac  Che  co¬ 
efflclenC  of  resclCuClon  remains  consCanc  over  Che  enclre  range  of 
veloclcles.  This  assumpclon  Is  noc  enclrely  Crue^,  buC  when  auComaClc 
weapons  are  under  conslderaclon  Che  range  of  impacC  veloclcles  Is 
narrow  and  Che  change  In  Che  coefflclenC  of  resClCuClon  Is  crlvlal. 

2.  PROCEDURE 


a.  Consider  Cwo  bodies  moving  along  Che  same  pach  In  space  aC 
consCanC  veloclcles  and  V2.  If  Che  veloclcles  are  such  ChaC  Che 
leading  body  will  noC  ouCrun  Che  oCher,  Chey  will  evenCually  collide 
wlch  a  resulclng  loss  of  energy.  Ic  Is  assumed  Chac  upon  colliding, 
a  force  Is  exerCed  on  Che  bodies  which  Is  analogous  Co  che  force 
developed  by  a  spring  wlCh  viscous  damping. 

b.  Figure  1  deplcCs  a  cyplcal  cwo  body  sysCem  aC  che  Clme  of 
InlClal  conCacc.  PqslClons  relaclve  Co  clme  are  given  In  recCangular 
coordinaces  aC  C  s  q,  by  xi  -  X2  -  0,  iti  -  Vi  and  X2  -  ±  V2. 


Figure  1 
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The  damped  spring  will  produce  a  force  affecting  the 
motion  of  the  bodies  only  when  the  relative  displacement  (X1-X2) 
is  positive*  The  equations  of  motion  are  as  follows: 


M  X,  +  R  k  (x,- Xi)  =  0 

■»-  R (  Xi-x.)  +  l<(Xi-x,)  =  o 

Initial  conditions:  =  X.z=  0  ,  X,  =  V,  ,  Xj^  = 

The  general  solution  of  this  equation  is  dependent  on  the  roots 
of  the  auxiliary  equation: 


This  equation  is  similar  to  that  obtained  for  one  degree  of 
freedom  where  the  viscous  damping  may  be  such  that  spring  may  be 

a.  under  damped  4KMqJ  R^(q+1) 

b*  critically  damped  4KMq  ;  R^(q+1) 

c*  over  damped  4KMq<R^(q+l) 


'^he  over  damped  and  critically  damped  cases  are  of  minor  im* 
portance  in  this  study  as  they  produce  a  coefficient  of  restitution 
of  zero.  They  are  considered  as  special  cases  in  Section  2  e 


The  under  damped  case,  which  produces  an  oscillatory  motion, 
is  of  interest  because  of  its  potential  to  define  a  coefficient  of 
restitution  in  the  range 

0<a  <  1  . 

Then, if,  4KMq^  R^(q+1) ,  the  general  solution  of  equation  (1) 


X,  =  A, +  5,t  + 

e  £  D,  sin  tot +  £/ coscotj 

(2) 

Xi  =  Az  +  bz'^ 

[  P^sin  u}t -f-E^cosut] 

wiere  O 

:  R 

(2a) 

\ 

and  ^ 

Z%M  ^ 

(2b) 

Applying 

the  initial  conditions  yields 

(3) 

is : 
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X;j  =  Vit  -r  r„t 

w(<b+l)  L 


-Pt  . 

e  sin  w 


*] 


To  obtain  the  separation  velocities  It  Is  necessary  to  solve  for 
the  time  i^en  the  relative  displacement  becomes  zero.  Hence,  subtract¬ 
ing  equations  (3)  and  equating  to  zero,  yields 


X 


XrXa 


UJ 


sin  u't  -  O 


(4) 


Let  t2  be  the  time  vdien  x  :  0,  then 
S  i  O  O 


(5a) 


and 


(5b) 


Differentiating  equations  (3)  and  substituting  equations  (5) 
yield  the  velocities  Immediately  after  separation. 


X 


<;  =  V. .  L 


J_ 

% 


-  e 


_  9*rr 


] 


Subtracting  to  obtain  the  final  relative  velocity 

X'-  x;-xi= 

As  the  Initial  relative  velocity  Is 

X  »  X,  -  X2 


(6) 


(7) 

(8) 


then  by  definition,  the  coefficient  of  restitution  may  be 
obtained  from  (7)  and  (8) 


0,  = 


=  e 


(9) 


Equation  9  defines  the  coefficient  of  restitution  as  a  function 
of  P  andu>  which  are  functions  of  M,  K,  and  R.  It  Is  Important  to  note 
that  the  coefficient  of  restitution  Is  completely  Independent  of  the  magnl 

tude  of  the  velocities. 
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Now  letting  -  _L  • 

L 

then  from  equations  (9),  (2)  and  (2b) 


/  L  ^  '77'  .  ‘77'/?  NTo+T 


Tnerefore 


R  =  /"  t  .  //-/w 


(10) 


(11) 


and 


U) 


■  »-  fWsrirJ 


(12) 


Equations  (11)  and  (12)  define  the  damping  factor  R  and  the 
frequency  of  oscillation  U>'as  a  function  of  K,  M,  q  and  a. 

c.  Maximum  Displacement.  When  two  bodies  collide  some  in¬ 
dentation  will  exist,  or  our  fictitious  spring  will  be  compressed. 

It  is  highly  desirable  to  have  a  spring  rate  which  is  realistic; there¬ 
fore  an  idea  as  to  the  possible  compression  of  the  spring  is  required. 


From  equation  (4)  we  have 

G  '  s  I  n  u)t! 


CJ 


Differentiating, 

X  =  V, -V, 


-ft 


CJ 


CJ  cos  CJ 


t- 


S\y)  CJ 


If  the  maxlmuin  displacement  occurs  at  t  =  t^^then 


.  tan  ct,  ^  _i_ 

P  In  b  -h  ! 


(4) 


(13) 


(14) 


or 


-l.r 

^  1}  In  b 


(15) 


U) 
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Substituting  tst]^  into  equation  (4)  yields: 

^max  *■  -Z3  [  (c^i-l)  /nb  J 

d.  Consistency.  The  derived  equation  1  thru  16  must  be  completely 
consistent  with  equation  A. 


If 

Ptt 


is  substituted  into  equation  (6),  the  equations  may  be  written  as 
follows : 


y'  *  V/*-  1/  +  "Vi  ~^3)  n 

(6a) 

yl  \/l  -  y  .  y.-VL  ,  CLf\/,-Vj 

Vz  -  Vi  * 

(6b) 

Substituting 

M=M,=  ^ 
into  equation  A  yields , 

V'  =  Vz-*-V,  Cl-aq,) 

(17) 

\/i  .  -(a-q,)  V^i. 

^  f) 

(18) 

When  the  terms  on  the  right  side  of  equations  (&a)  and  (6b)  are 
placed  over  a  common  denominator  and  the  terms  collected,  the  results 
are  identical  to  those  of  equation  A.  It  is  evident  that  the  deriva¬ 
tions  are  mathematically  consistent. 

e.  Special  Cases 

(1)  One  Degree  of  Freedom.  The  equations  1  through  16  may  be 
readily  utilized  to  represent  a  system  in  one  degree  of  freedom.  If 
the  mass  of  the  second  body  is  assumed  to  be  infinite  and  motionless, 
q  will  equal  00  and  V2  will  equal  0  ; 
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therefore 


j-ti-  and  j 

%  -jr  -! 

(2)  Critical  Damping.  A  special  case  of  critical  damping  occurs 

when 

4  K  Mq  .  (q+1) 


The  general  solution  of  equation  (1)  is  then 

-pt  -pt 

Xj  s  Aj  +  Bj^t  +  Di  e  +  Ejte 

-pt  -Pt 

X2  s  A2  +  B2t  +  D2  e  +  E2te  ' 

applying  initial  conditions 

XI  •  (q  V2  +  Vi)t  +  (Vi  -  V2)qte"^*^ 
q+1 

X2  «  (q  V2  +  Vi)t  -  (Vi  -  V2)t  e"P*^ 

the  relative  displacement  and  velocity  are  then 
X  s  X1-X2  «  te  ^'^l  “ 

X  ■  s  (Vj-V2)(l-pt)e 


(19) 
(19a) 

(20) 

(20a) 


(21) 

(21a) 


Figure  2 


-8- 


REPORT 

SA-TR20-2811 

The  relative  displacement  and  velocity  as  a  functions  of  time  are 
shown  in  Figure  2.  Since  this  is  a  case  of  critical,  damping,  the  rela¬ 
tive  displacements  and  velocities  will  be  zero  when  t  s® .  It  can  be 
seen  that  the  sign  of  the  displacement  is  never  changed;  consequently 
the  coefficient  of  restitution  is  zero. 

(3;  Over  Damping.  An  over  damped' condition  occurs  when 
r2  (q+l)>  4KMq 

The  general  solution  of  equation  (1)  is  then 

X,"  A,  +  5,t+ 

where  -Pj^  and  -P2  are  the  negative  real  roots  of 

The  over  damped  case,  like  the  critically  damped  case,  produces  a 
coefficient  of  restitution  equal  to  zero.  The  difference  in  the  two  cases 
is  that  the  velocity  of  the  overdamped  system  becomes  zero  at  some  finite 
time  and  the  displacement  does  not  approach  zero. 


It  may  readily  be  seen  that  as  the  relative  displacement  does  not 
change  sign,  the  two  masses  are  lumped  into  one  and  the  original  masses 
are  permanently  deformed. 

3.  DISCUSSION 


a.  The  derived  equations  may  now  be  applied  to  a  dynamic  system 
to  evaluate  the  rate  and  damping  of  a  fictitious  spring  and  the  effects 
of  the  coefficient  of  restitution. 

These  equations  indicate  that  some  displacement  must  be  tol¬ 
erated.  In  the  case  of  solid  collisions,  a  very  small  displacement  is 
permissible*  In  the  case  of  buffers,  the  displacement  obtained  math¬ 
ematically  will  be  lower  than  that  observed  because  the  derived  equa¬ 
tions  do  not  take  into  account  the  rate  of  the  buffer  spring.  In 
either  case,  the  small  displacements  during  collision  are  not  signifi¬ 
cant  if  the  velocities  and  time  of  contact  are  correct. 
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b.  M-73  Machine  Gun  Analysis,  In  the  M-73  machine  gun  two  ma¬ 
jor  collision  systems  exist,  in  which  the  coefficient  of  restitution 
is  important.  The  characteristic  type  of  each  is  described  as  follows: 

(1)  Solid  Stop.  A  collision  between  barrel  extension  and 
receiver. 


(2)  Buffer.  A  collision  between  barrel  extension  and  buffer 
spring. 

c.  Time  displacement  curves  were  plotted  for  solid  stop  and  buf¬ 
fer  collisions  in  the  M-73  machine  gun  on  a  rigid  mount.  These  curves 
reflect  the  following  conditions: 


Solid  Stop  Collision 


q  .  00 

M  «  9.5/386. 


#sec^ 

in 


Buffer  Collision 

q  =  00 

M  -  9.5/386 


V  ■  60  in/sec 


V  s  55  in/sec 


V  »“5  in/sec 


V  r-30  in/sec 


t2  =  0.002  sec 

with  ®  ■  1 

b 


t2  »  0.00254  sec 


Applying  equation  (10) 

b  -  12  b  «  11 

6 

In  b  -  2.48491  In  b  =  0.60014 

From  equations  (5)  and  (12) 

K  =  98,616.25  K  -  39,011.35 

Substituting  into  equation  (12) 
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z  1.570.795 

Then  from  equation  (11) 
R  a  In  b 


=  1,236.846  . 


4  /CM 

Solving  for  from  equation  (16) 

-I 


^ Wdx  (jj  ®  sin  ^'ta.n 


tr 

In  b 


} 


ti  S  1.20426 
Substituting  into  equation  (15) 

CJ'tf  ^ 

//>  O 

tj^  «  0.57393  m  sec. 

and 

sin  Urt^  =  0.78431 
Then  from  equation  (10) 

In  b-  e-^ 

P  =  1.24245x1o3 
^ti  =  0.71308 
e'P*^!  .  0.49013 


t^  .  5.18294  . 


t]^  s  1.11589  m  sec. 


sin  ur't. 


0.98189 


p  :  0.23864x103 
^t^  s  0.26630 
e’P*^!  >  0.76621  . 


Substituting  back  Into  equation  (16) 
-Pt, 


X 


max 


-3 

X  max  •  14.7x10  in  x  max  s  33.45x10  in 
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d.  Briefly  stated,  the  summarized  results  of  the  preceding  calcu¬ 
lations  are  as  follows: 

(1)  Solid  Stop  Collision^  With  K  -  98,600  #/in  and  R-61.1,  a 
coefficient  of  restitution  of  0.08333  will  be  maintained 
if  a  maximum  displacement  of  14.7  x  10"^  inches  is  accept¬ 
able. 

(2)  Buffer  Collision.  With  Ks39,000  and  R-11.73  a  coefficient 
of  restitution  of  0.54545  will  be  maintained  if  a  maximum 
displacement  of  33.45  x  10*^  inches  is  acceptable. 

3  4 

e.  Analog  Computer  Application. *^The  analog  computer  program  for 
the  M-73  machine  gun  is  such  that  actual  time  is  slowed  down  by  a 
factor  of  200.  Hence,  the  forces  simulating  losses  due  to  the  coef¬ 
ficient  of  restitution  will  be  of  very  short  duration. 

To  check  the  utility  of  the  derived  relationships,  consider  a  con¬ 
dition  during  which  the  barrel  extension,  moving  with  free  velocity, 
strikes  the  receiver  and  bounces  back.  This  condition  may  be  satisfied 
by  the  following  differential  equation: 

Mx  Fg  a  0 

with  initial  conditions 


taO  x«0.600  X  *  -60  in/sec 

where  Fg«0  when  x  ^  0 

and  Fs  ■  98,600x  +  61.1  x  when  x<0 
Substituting  the  values  of  mass  and  dividing 

X  +  2,485. 162^  +  4,010,425.26x  a  0 
(1)  Computer  Equations. 

i-P  T  -  zoo  t 

X  =  2  X 
X  =  400 X 
^  =  8  0^000  X 

^  ^  /2.425'Sl  X  +  lOO-ZCS.  “  0 

uK.n  T'O;  X  =  30.  OO  vo/fs  ^  X  =  iS.OO^^olts 
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f.  Patch  Dtagram.  The  patch  diagram  illustrated  in  Figure  3  is 
used  to  program  the  analog  computer  for  evaluation  of  the  collision 
systems  in  the  M-73  machine  gun. 


This  computer  program  has  been  tested  repeatedly  with  highly  satis¬ 
factory  results.  The  coefficient  of  restitution  was  evaluated  for  a 
range  of  velocity  factors  of  4  to  100  volts.  Slight  deviations  in  the 
coefficient  of  restitution  were  noted  when  the  velocity  was  below  3 
volts . 
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